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Abstract 

In order to understand the parameters of the standard model of 
electroweak and strong interactions (coupUng constants, masses, mix- 
ing angles) one needs to embed the standard model into some larger 
theory that accounts for the observed values. This means some addi- 
tional sector is needed that fixes and stabilizes the values of the funda- 
mental constants of nature. In these lecture notes we describe in non- 
technical terms how such a sector can be constructed. Our additional 
sector is based on rapidly fluctuating scalar fields that, although com- 
pletely deterministic, evolve in the strongest possible chaotic way and 
exhibit complex behaviour. These chaotic fields generate potentials 
for moduli fields, which ultimately fix the fundamental parameters. 
The chaotic dynamics can be physically interpreted in terms of vac- 
uum fluctuations. These vacuum fluctuations are different from those 
of QED and QCD but coupled with the same moduli flelds as QED 
and QCD are. The vacuum energy generated by the chaotic flelds 
underlies the currently observed dark energy of the universe. Our 
theory correctly predicts the numerical values of the electroweak and 
strong coupling constants using a simple principle, the minimization 
of vacuum energy. Implementing some additional discrete symmetry 
assumptions one also obtains predictions for fermion masses, as well 
as a Higgs mass prediction of 154 GeV. 



1 Introduction 



String theories predict an enourmous amount of possible vacua after com- 
pactification, of the order 10^^°. In each of these vacua the cosmological 
constant as well as the fundamental constants of nature can have different 
values. One is led to the so-called 'landscape' picture p. The landscape 
in a sense represents the set of all possible blueprints of the universe. To 
select the right vacuum, i.e. the one we observe around us right now, often 
an anthropic point of view is chosen. But is an anthropic principle really the 
last word and the ultimate answer to all our questions? 

A natural idea to avoid an anthropic selection principle would be that 
the universe will not be left alone with its choice of 10^^° vacua but that it 
gets some help. This help should have the form of an additional sector, a 
theory as yet not included in the ordinary standard model, neither in ordinary 
string theories. The additional sector should yield some general principle to 
select, fix and stabilize the standard model parameters in the way we do 
observe them. In fact, it should do much more: It should choose the right 
gauge groups, the correct amount of supersymmetry, the correct number of 
flavour families, it should create an excess of matter over anti-matter, and 
so on. The additional sector should fix the most relevant information about 
the future universe at a very early stage, similar as a DNA string fixes the 
most important information of a human being already when the first cells 
are formed. 

We are still far away from such a theory. But some numerical observa- 
tions have recently been made [21 [Sj H] that seem to give us a hint how this 
additional sector could look like. Our aim here is to explain the relevant 
concepts in a non-technical way. 

In practice standard model parameters can be thought of as being fixed 
by so-called moduli fields. A varying standard model parameter (e.g. the 
fine structure constant) can be essentially regarded as (a simple function of) 
such a moduli field. These moduli fields evolve to minima of their potentials. 
So if we know the correct moduli potentials describing the world around us, 
we also know the correct standard model parameters. So what could be 
a theory to construct these moduli potentials? In principle the potentials 
should follow from the embedding theory (e.g. M theory + compactification 
+ supersymmetry breaking), but little is known in practice due to the enour- 
mous complexity inherent in the above theories. But as with any unknown 
theory we can be guided by first trying to find an empirical theory that does 
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the correct thing, i.e. reproduces the observed value of the fine structure 
constant and other fundamental constants, and then later try to embed it 
into a greater context. The interesting thing is that such an empirical theory 
is possible [21 |3l |3] : There is a class of highly nonlinear chaotic dynamical 
system that seem to reproduce the 'correct' standard model parameters by 
a simple selection mechanism, the minimization of vacuum energy. 

Physically the above chaotic dynamics can be regarded as describing 
rapidly fiuctuating scalar fields associated with vacuum fiuctuations. These 
are vacuum fiuctuations different from those of QED or QCD, so what would 
be the physical embedding? The most natural embedding is to associate 
the above chaotic vacuum fiuctuations with the currently observed dark en- 
ergy in the universe [U [S]. Since nobody really knows what dark energy 
is there is a lot of freedom in the dark energy sector and enough 'space' 
to embed new things! The chaotic fields, living in the dark energy sector, 
generate effective potentials for moduli fields — just the same moduli fields 
that are responsible for the fundamental constants of the standard model 
of electroweak and strong interactions. The moduli fields then move to the 
minima of the potentials generated by the chaotic fields, and fix the funda- 
mental constants of nature. The chaotic sector appears to provide a possible 
answer to the question why we do observe certain numerical values of stan- 
dard model parameters (such as the fine structure constant or the strong 
coupling at I^-mass scale) in nature, others not. It can be used to avoid 
anthropic considerations for fundamental constants. Moreover, it generates 
a small cosmological constant in a rather natural way [1] . 

These lecture notes are organized as follows: In section 2 we provide some 
background information on moduli fields, variable fine structure constants 
and related topics. In section 3 recall the method of stochastic quantization 
introduced by Parisi and Wu [6, ,7j. We then introduce the chaotic fields in 
section 4, and show how they can generate potentials for moduli fields in 
section 5. Finally, in section 6 we provide numerical evidence that there are 
various local minima of the potentials that do reproduce known standard 
model parameters with high precision. We obtain excellent agreement with 
measured values seen in collider experiments. We will also make a few pre- 
dictions for unknown parameters of the standard model, such as the Higgs 
mass, based on the chaotic field dynamics and some additional symmetry 
assumptions. 
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2 Moduli fields and variable fine structure 
constant 

The action of the standard model is well known [8], we don't have to work 
that out here. The term of the action involving the electromagnetic field 
tensor F^^, is given by 

Sp = J d'x^ {-\f,,f>^^^ . (1) 

In models with a variable fine structure constant a(t) [9l [TOl [TTl [121 [13] the 
above action is modified as follows: 

Sf = J d'xy^ ^-^-Bf{x/M*)F,,F^''^ . (2) 

Here x{t) is a homogeneous scalar field, henceforth called the moduli field (a 
related field in string theory is the dilaton field). M* is a mass of the order of 
magnitude of the reduced Planck mass, i^i? is a function that is in principle 
determined by the embedding theory (e.g. string theory or M theory). In 
Bekenstein's model [I2|, 

S^(X/M*) = e-2(>^-*)/*^*, (3) 

but other choices are possible as well. Changes in the field x{t) imply changes 
in the fine structure constant a{t). The relation is 

^ = I (4) 

a BF{x{t)/M*) ^ ^ 

where Bp{x/M*) = 1. Here d is the stationary value of a{t), and x the 
stationary value of xit)- 

The interpretation of the above generalized action ([2D is quite obvious: 
Normally, the electromagnetic field tensor would contribute with the term 
— ^F^yF^^ in the action. This relation is now only satisfied in the stationary 
case. Otherwise the strength of this term is given by a prefactor Bp that 
depends on the value of the field xif) at a given time t. The fine structure 
constant a{t) is a simple function of the moduli field x(t) according to eq. 

For small displacements of the field x(^) from its stationary value we may 
write 

Bp{x{t)/M*) ^ Bp{x/M*) + ^^^^tlAB'p{x/M*), (5) 
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where B'p is the derivative of the function Bp. We thus obtain in leading 
order from eq. (jlj) 



c^{t) = a[l-B'p{x/M*)^^^^Y (6) 

In Bekenstein's model, Bp{x/M*) = —2. Clearly, the fine structure con- 
stant approaches its stationary value d if the moduli field x approaches its 
stationary value x- Using Bp{x/M*) = —2 we may also write eq. ([6]) as 

hr^^^x(t)-x (7) 

/ a 

which shows that the fine structure constant and the moduli field are basically 
the same. 

To introduce a dynamics for the moduli field Xy need to know its 
potential V{x)- ^ Robert son- Walker metric the dynamics is then given by 

where H is the Hubble parameter and is the coupling of the field x to 
matter (in particular dark matter), is the matter density of the universe 
at a given time. In typical models studied in the literature [9], the coupling 
C,m is very small. 

Whereas in general the potential V is unknown, near to the equilibrium 
point X we may expand it as 

V{x) = ^m'^{x- xf + const (9) 

so that 

-^=m''ix-x)- (10) 

In [9] m is chosen as an extremely small mass parameter, of the order of the 
current value Hq of the Hubble parameter: 

mr^ Ho (11) 

That is to say, one considers an ultralight scalar field x with a mass of 
the order 10~^^ eV. There is motivation for the existence of such ultralight 
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fields from extended supergravity tlieories [2]. In terms of the fine structure 
constant a{t), one obtains near tlie stationary point 



V{x) = -m'M- 



1 



(a(t) -af 



+ const 



(12) 



(S^(X/M*))2 



By construction, tlie energy density m?'M*'^ associated witli tliis potential is 
of the same order of magnitude as the dark energy density of the universe at 
the present time. 

Of course, in a similar way one can introduce further moduli fields cor- 
responding to other standard model coupling constants, masses and mixing 
parameters as well. The fine structure constant was just one example. In fact, 
for each relevant standard model parameter there should be a corresponding 
moduli field. So we need about 20 such fields. 

3 Stochastic quantization 

Let us now proceed to a 2nd-quantized theory. An elegant method to do 
2nd quantization is via the so-called stochastic quantization method. In the 
Parisi-Wu approach of stochastic quantization one considers a stochastic dif- 
ferential equation evolving in a fictitious time variable s, the drift term being 
given by the classical field equation [6l [7] . Quantum mechanical expectations 
correspond to expectations with respect to the generated stochastic processes 
in the limit s ^ oo. The fictitious time s is different from the physical time 
t, it is just a helpful fifth coordinate to do 2nd quantization. Neglecting 
spatial gradients the field under consideration is a function of physical time 
t and fictitious time s. The stochastically quantized equation of motion of a 
homogeneous scalar field (p in Robertson- Walker metric is 



where H is the Hubble parameter, V is the potential under consideration 
and L{s,t) is Gaussian white noise, 5-correlated both in s and t. For e.g. a 
numerical simulation we may discretize eq. f|T3|) using 




(13) 



s 



t = i6, 



(14) 
(15) 
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where n and i are integers and r is a fictitious time lattice constant, 5 is a 
physical time lattice constant. The continuum limit requires r ^ 0, 5 — 0, 
but since quantum field theory is expected to break down at the Planck scale, 
it can make physical sense to keep a small lattice constant of that size as an 
effective cutoff: 

1,1 

r ~ — 5-, b ~ . (16) 

mj^i mpi 

Eq. f[T^ yields the discrete dynamics 

= ^M^'-2vl. + Vn') + 3§{vl-^l') + V\^l,)+nozse. (17) 
This can be written as the following recurrence relation for the field y?^ 

= (1-a) jy^t. + ^\/'(^D|+3— (y.j,-(^r') + 2(¥'n^'+¥'n"')+™^^e, 

(18) 

where a dimensionless coupling constant a is introduced as 

a := |. (19) 

This coupling a is a free parameter of our stochastically quantized theory. A 
priori it can take on any value, reminescent of a moduli field. 

We may also introduce a dimensionless field variable $^ by writing = 
^nPmax, whcrc Pmax IS some (so far) arbitrary energy scale. The above scalar 
field dynamics is equivalent to a spatially extended dynamical system (a 
coupled map lattice) of the form 

K+i = (1 - a)n^'J + ^HSai^^^ - $-1) + ^($^,+1 + $-^) + r ■ no^se, (20) 
where the local map T is given by 

T($) = $ + ^r^V'iPma.^). (21) 



Here the prime means 

d 1 d 



(22) 



= (1 - a)T(<|.l) + + $^1) + r ■ noise (23) 



dlf Pmax <9$ 

A symmetric dynamics of the form 

a 
2 



is obtained if H5 << 1, equivalent to 



5 « H-'. (24) 

This approximation is valid if the universe is much older than the physical 
time lattice constant 5. In this case the term proportional to H in eq. fl2Ul) 
can be neglected. The local map T depends on the potential V under consid- 
eration. Since we restrict ourselves to real scalar fields T is a 1-dimensional 
map. 

Let us summarize the main result of this section: Iterating a recurrence 
relation of the form 0231) is equivalent to considering a stochastically quantized 
scalar field. The relation between the map T and the potential V is 

V'{^) = ^ l-^+Pma.T (^]]. (25) 
T I \Pmax I J 



Integration yields 

= l^|_i<^2^p^^^| d^T^y-^-^^ const, (26) 

which in terms of the dimensionless field $ this can be written as 

V^(^) = + I dm^)^ + const. (27) 



4 Introducing chaotic fields 

We now come to the crucial point, namely what type of dynamics is generated 
by fl25]) for various types of scalar field theories. Take the example of an 
ultralight moduli field x- For these types of fields Pmax ~ ^* ~ is large, 
whereas V'ix) ~ "^ix ~ X) is extremely small. This means the mapping T 
in eq. (!2T]) is extremely close to an identity and the field moves very slowly. 
Effectively this means that our discrete dynamics (1231) approximates very 
well a smooth continuum evolution of the field x- The field x smoothly 
approaches a minimum of the potential and stays there, apart from some 
small fluctuations induced by the noise term. 

Once again, for a moduli field x the energy scale associated with the 
field variable is very large, of the order of the Planck mass M* ~ nripi, 
whereas the potential V contains an extremely small mass parameter m, of 
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the order of the present Hubble constant Hq. We could ask whether for 
symmmetry reasons maybe another scalar field (p exists that has just the 
opposite properties, i.e. the energy scale of that field variable ip is extremely 
small (of the order Hq) but its potential V contains a mass term of order 
M*? 

Indeed, such a field has been studied in Due to the fact that the 
forcing is now very strong, this field exhibits strongly chaotic behaviour if 
it is stochastically quantized. This can immediately be seen from eq. (!2T]) : 
With Pmax being small, the map T is now far away from the identity. Chaotic 
behaviour is possible. The small noise term in eq. (!23|) can actually be 
neglected in the chaotic case. 

We may go one step further in our symmetry considerations. The moduli 
field is as slow and regular as a field can be. What would be the other 
extreme? What would be a field that is as rapidly fluctuating and irregular 
as possible? In other words, which scalar field dynamics would create the 
strongest possible chaotic behaviour? 

The above question has been solved in [TSl [TS] and is well understood 
from a nonlinear dynamics point of view. One knows that the maps with 
strongest chaotic properties (being smooth and deterministic at the same 
time) are those conjugated to a Bernoulli shift [17] (a shift of integer sym- 
bols in suitable coordinates). Among those, certain deterministic maps are 
even more 'random' than others, in the sense of having least higher-order 
correlations: These are the so-called Tchebyscheff maps T/v of A^-th order, 
defined as 

T2($) = 2$2 - 1 (28) 
T3($) = 4$^ - 3$ (29) 

(30) 

T^($) = cos(Ararccos$) (31) 

They can arise out of the above dynamics (|23|) for suitable potentials V. 

The most important scalar field potential in particle physics is of course 
a double-well potential. So let us consider the distinguished example of a y^^- 
theory generating strongest possible chaotic behaviour. Take the potential 

V.M = 1^ _ + const, (32) 

^ t Pmax ) 
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or, in terms of the dimensionless field $, 

V^M = ^^pLA^' - + const. (33) 

Applying our formalism of the previous section, we end up with the following 
local map: 

$n+i = T_3($n) = -4$^ + 3$„ (34) 

T_3 is the (negative) third-order Tchebyscheff map. It is conjugated to a 
Bernoulli shift of 3 symbols, and generates the strongest possible stochastic 
behaviour possible for a smooth low-dimensional deterministic dynamical 
system. 

Apparently, starting from the potential ( !32l) we obtain by second quan- 
tization a field ip that rapidly fluctuates in fictitious time on some finite 
interval, provided that initially ipo G [—pmax,Pmax]- The small noise term in 
eq. (123!) can be neglected as compared to the deterministic chaotic fluctua- 
tions of the field. We physically interprete these rapid changes of the field (p 
as representing vacuum fluctuations. Of course these are vacuum fluctuations 
different from those of QED or QCD. Since the expectation of the vacuum 
energy associated with the chaotic field is (V^s^ip)) ~ Pmax/'^ ~ -^o"^p^ ^^^^ 
a chaotic field yields the correct order of magnitude of vacuum energy density 
in order to account for a small cosmological constant. Hence we assume that 
the chaotic fields underly dark energy. 

The above example of a chaotic yj'^-theory can be generalized. There are 
various discrete degrees of freedom to introduce a deterministic chaotic field 
dynamics that is as random as possible. Consider a 1-dimensional lattice, 
the lattice sites are labelled by integers i. At each lattice site i we have 
a dimensionless field variable which evolves in discrete time n. In the 
uncoupled case, the dynamics is given by 

= ±T^m, (35) 

where iT^v is either the positive of the negative Tchebyscheff map. The 
initial value $o is chosen on the interval [—1,1], the iterates $^ then stay 
in this interval, but evolve in a deterministic chaotic way. The dynamics is 
conjugated to a Bernoulli shift of symbols, which means that in suitable 
coordinates the iteration process is like shifting symbols in a random symbol 
sequence. There are further discrete degrees of freedom to do the coupling 
to the nearest neighbours. Instead of coupling to the variables ^^^^ at the 
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neighboured sites as in eq. (1251) . we could also use the updated variables 
±Tjsf{^'^^^). Since n describes fictitious time, it's really not clear which choice 
is the correct one. For sure both degrees of freedom exist. All these discrete 
degrees of freedom can be written in a compact form as follows: 

K+i = (1 - a)Tr,{K) + s^{TUK-') + TUK^')), (36) 

where s is a sign variable taking on the values ±1. The choice s = +1 
is called 'diffusive coupling', but for symmetry reasons it also makes sense 
sense to study the choice s = — 1, which we call 'anti-diffusive coupling'. The 
integer b distinguishes between the forward and backward coupling form, 
6 = 1 corresponds to forward coupling (T^($) := T/v($)), 6 = to backward 
coupling (T^($) := $). We consider random initial conditions, periodic 
boundary conditions and large lattices of size imax- 

One can easily check that for odd the choice of s is irrelevant (since 
odd Tchebyscheff maps satisfy T/v(— $) = — T/v($)), whereas for even the 
sign of s is relevant and a different dynamics arises. Hence, restricting our- 
selves to A^ = 2 and A^ = 3, in total 6 relevant chaotic scalar field theories 
arise, characterized by (A^, 6, s) = (2, 1, +1), (2, 0, +1), (2, 1, -1), (2, 0, -1) 
and {N,b) = (3,1), (3,0). For easier notation, in the following we will la- 
bel these chaotic field theories as 2A, 2B, 2A^ , 2B~ , 3A, 3B, respectively. 

The important thing to remember from this section is that eq. fl36|) just 
describes a degenerated stochastically quantized scalar field with strongest 
possible chaotic properties. For somebody from the nonlinear dynamics and 
complexity community, the dynamics (l36ll represents the most natural thing 
in the world: It's just a coupled map lattice, a standard example of a spatio- 
temporal dynamical system [18] . For somebody from the elementary particle 
physics or string theory community, however, eq. fl36p may look somewhat 
unusual and unfamiliar at first sight. But it's sometimes worth to learn new 
things! 

The dynamics fl36|) is a dissipative deterministic dynamical system. It 
exhibits chaotic behaviour and produces information, measured by a positive 
KS entropy [T7]. Dissipative deterministic systems as a model of quantum 
gravity have also been suggested in [I9j . 
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5 Generating potentials for moduli fields 



Clearly, if the coupling a in the above chaotic field dynamics (136|) is chosen 
as a = 0, then there are no correlations between neighboured lattice sites: 



Here the notation (■ ■ ■) denotes the expectation value, which can be nu- 
merically calculated by iterating the maps for random initial conditions and 
doing a time average. Note that the index i denotes physical time (in units 
of the lattice constant 5). If we physically interpret the chaotic fiuctuations 
of the field $^ as some sort of vacuum fiuctuations underlying the dark en- 
ergy of the universe (see [1] for more details), then the above condition (1571) 
looks physically very reasonable: We want subsequent vacuum fiuctuations 
to be uncorrelated in physical time, because otherwise they wouldn't really 
describe spontaneous fiuctuations. 

We may, however, insist on an interacting theory, i.e. a 7^ 0. Are there still 
some distinguished couplings a* 7^ where we can keep the above condition 
of a vanishing correlation function in physical time? Of course such a state 
would again be distinguished as being as random as possible, where this 
concept is now extended from fictitious time to physical time. There are also 
a couple of other reasons why one wants a vanishing correlation function, see 
[3] for more details. 

In fact, numerical investigations show that the above states with vanish- 
ing correlation exist and distinguish certain coupling constants a* . These 
are numerically observed to coincide with known standard model coupling 
strengths (see next section). In general one can show [21 [3] that the quantity 



can be physically interpreted as the interaction energy of the chaotic field 
theory under consideration (in suitable units). Hence states of strongest ran- 
dom properties, described by a vanishing correlation function, have vanishing 
interaction energy. 

Besides the interaction energy, there is also another relevant vacuum en- 
ergy associated with the chaotic fields. This is the self energy V{a), given 




) = 



(37) 




(38) 



by 




(39) 
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respectively 

r(3)(a) = ^($2) _ ($4) (AT = 3). (40) 

For a derivation, see The self energy describes the vacuum energy as- 
sociated with the potentials that generate the chaotic dynamics in fictitious 
time (the additive constant is fixed by some symmetry considerations). 

In the following, we will use both the interaction energy and the self 
energy to generate suitable potentials for moduli fields. Recall that classically 
the moduli field x obeys 

X + 3i7x + ^'(x)=0, (41) 

where we neglected possible interactions with dark matter. In the vicinity of 
the stationary state Xi oiis has V\x) = "^^(x ~ x) and the moduli field is 
essentially the same as the standard model coupling constant a. According 
to eq. (171) we have 

x{t) = ^a{t) + x-lM*. (42) 
Putting eq. fjl2|) into fHTl) . we obtain an equation for a, 

a + 3Ha + V'{a) = 0, (43) 
where locally the potential is given by 

V{a) = ^m'{a-af. (44) 

Eq. is just the equation of a damped harmonic oscillator, provided a is 
in the vicinity of a. 

The crucial point is now that we assume that the same moduli fields 
are responsible for the coupling constants of the standard model of elec- 
troweak and strong interactions and those of the dark sector described by 
the chaotic fields. This means a = a. The chaotic field wants to find a state 
of strongest possible random properties described by a vanishing correlation 
function ($^$^"'"^). Hence the moduli field x given by eq. fH2|) with a = a 
needs to adjust its value. This can be achieved by choosing in eq. fH3l) the 
formal potentia0 

V{a) = -m^ r dol^^'W^^). (45) 
h 



maps. 



^Our sign conventions relate to moduli potentials generated by positive Tchebyscheff 
,r)s. 
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Figure 1: Basic idea underlying this paper. A priori all kinds of standard 
model couplings a = a (values of moduli fields) are possible. These are then 
driven into stable zeros of the interaction energy W{a) of the chaotic fields. 
In the above schematic picture, the stable zeros are denoted by Oi and 02. 

By construction, local minima of this potential are stable stationary points 
corresponding to stable (attracting) zeros of the forcing V'(a) = — m^($^$^"^^) 
(see Fig. 1). If the damping H in eq. 0431) is large as compared to m, and if 
the initial displacement is not too large, then the stationary state a is rapidly 
approached. 

For other types of moduli fields, we may choose in eq. ( I43l) the self energy 
potential V{a) = rn^V^^^a) as given in eq. ( l39l) or ( HOl) . These types of 
moduli fields then evolve to local minima of the self energy. 

It is well known that standard model interaction strengths actually de- 
pend on the relevant energy scale E. We have the running electroweak and 
strong coupling constants. For example, the fine structure constant aei{E) 
slightly increases with E, and the strong coupling rapidly decreases with 
E. 

What should we now take for the energy (or temperature) E of the moduli 
fields near to their stationary state? In other words, if a standard model 
coupling a is fixed as a minimum of the potential V{a), at which energy scale 
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E is this standard model coupling fixed? A priori this is unknown. However, 
there is extensive numerical evidence ^ ^ that the distinguished couplings 
a = a* which correspond to minima of the potentials numerically coincide 
with running standard model coupling constants a{E) with an energy E 
being given by 



Here N is the index of the Tchebyscheff map considered, and mB^rnf^,mj^ 
denote the masses of a boson B and a fermion fi and anti-fermion /2 — not 
just some exotic bosons and fermions but precisely those that we know from 
the standard model. One typically observes particle combinations B,fi,f2 
that describe possible interaction states in the standard model, for example 
a decay of the form 5 — >■ /i + /2 or a reaction /i + /2 — >■ B. Detailed 
evidence will be given in the next section. Formula (146!) formally reminds us 
of the energy levels Epf = ^ quantum mechanical harmonic oscillator, 

with low-energy levels {N = 2,3) given by the masses of the standard model 
particles. 

6 Numerical results 

We now present our numerical results (much more details can be found in 
[2], [3]). In particular, we will show that stable zeros of the correlation func- 
tions (states of strongest random properties for chaotic fields, which are sta- 
tionary states for moduli fields) reproduce known standard model coupling 
constants and masses with high precision. Our results are obtained by long- 
term iteration of the chaotic dynamics. One numerically calculates, for a 
given coupling a, the interaction energy as the time average 



where the $^ evolve according to (136!) with random initial conditions. In 
practice, we used a finite lattice of size J = 10000 with periodic boundary 
conditions, and iteration numbers corresponding to several weeks of CPU 
time. Everybody is welcome to reproduce and verify the numerical results 
described below — the recurrence relation (|36l) can be very easily installed on 
any computer. What one observes is that the stable zeros of W{a) coincide 



E = -N{mB + rrif^ + mjj. 



(46) 



1 1 1 M J 

Wia) = = - lim lim — - V V <l>t$i+\ 



(47) 
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Figure 2: Interaction energy of the 3A dynamics in the small-couphng region. 



with known standard model interaction strengths, thus indicating the phys- 
ical relevance of the theory presented presented in sections 1-5. The chaotic 
fields appear to determine the moduli potentials in precisely the way we 
need them to be for a physically realistic vacuum. They can 'help' the uni- 
verse to find the 'right' vacuum out of an incredibly large number of choices. 
Anthropic considerations can be avoided in this context. 



6.1 The 3 A dynamics — electric interaction strengths 
of electrons and d-quarks 

Fig. 2 shows the interaction energy W{a) = |($^$^"^^) of the chaotic 3A 
dynamics in the small-coupling region. We observe two stable zeros of the 
interaction energy in the low-coupling region: 

ap^) = 0.0008164(8) 
4^^) = 0.0073038(17) 

Previously these stable zeros were denoted by a* or a. A stable zero satisfies 
V'{a*) = and V"{a*) > 0, with V given by eq. 
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Remarkably, the zero 02"^^ appears to numerically coincide with the fine 
structure constant aei ~ 1/137. To construct a suitable physical interpre- 
tation in the sense of eq. fH^ . let us choose i? = 7, /i = e~,/2 = e"*". 
The relevant energy scale underlying this moduli field is then given hy E = 
(3/2)(m^ + 2me) = 3me, according to eq. (HGll . Hence our standard model 
interpretation of the stationary moduli state described by 02'^'' suggests the 
numerical identity 

4'^^ = aeKSme). (48) 

For a precise numerical comparison let us estimate the running electromag- 
netic coupling at this energy scale. We may use the Ist-order QED formula 



a,i{E) = a,i{Q) 1 + '''' ' 



^Y^fi J dx x{l — x) log 1^1 H 2^(1 ~ ^) j I • 



(49) 

The sum is over all charged elementary particles, rrii denotes their (free) 
masses, and fi are charge factors given by 1 for e, /i, r-leptons, | for u, c, t- 
quarks and | for rf, s, 6-quarks. Using this formula, we get aez(3me) = 

0.007303, to be compared with a^^^ = 0.0073038(17). There is excellent 
agreement. Inverting the argument, the above zero of the interaction energy 
of the chaotic field can be used to predict the numerical value of the fine 
structure constant from first principles. 

Next, we notice that the other zero a^^^^ has approximately the value 
^aei- This could mean that the chaotic 3A field also has a mode that provides 
evidence for electrically interacting (i-quarks. Our interpretation is 

af^^ = aii{3md) = ^-aeiiSnid), (50) 

where af; = ^aei denotes the electromagnetic interaction strength of d- 
quarks. In the harmonic oscillator interpretation (l46l) . we may choose 5 = 7, 
/i = d, f2 = d. Formula (H9|) . as an estimate, yields for ^ 9 MeV the value 
aeiiSnid) = 0.007349, which coincides very well with 9af^^ = 0.007348(7). 
The value 90^'^'^'' actually translates to the energy scale Ed = (26.0 ± 6.4) 
MeV. This yields nid = ^Ed = (8.7 ± 2.1) MeV, which coincides with es- 
timates of the MS current quark mass of the d quark at the proton mass 
renormalization scale. 
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Figure 3: Interaction energy of the 3B dynamics in the low-couphng region. 

6.2 The 3B dynamics — weak interaction strengths of 
neutrinos and u-quarks 

The interaction energy W{a) of the SB field is plotted in Fig. 3. In the low- 
couphng region a e [0, 0.018] we observe the following stable zeros of W{a): 



ap^) = 0.0018012(4) 
: 0.017550(1) 



"2 



If our approach is consistent, we should be able to find an interpretation of 
ap^^ and a^^^ in terms of moduli fields fixing the standard model couphng 
strengths of w-quarks and neutrinos. 

Let us start with a2 ■ For left-handed neutrinos, the weak coupling due 
to the exchange of Z°-bosons is given by 



1 



(51) 



4 sin^ 6w cos^ 

Here 9^^ is the weak mixing angle. In the following we will treat sin^ 9w 
as an effective constant, and regard cte/ as the running electromagnetic cou- 
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pling. Other renormalization schemes are also possible, but yield only mi- 
nor numerical differences. Experimentally, the effective weak mixing angle 
is measured as sin^ 9w = sf ~ 0.2315 |20]. Assuming that in addition to 
the left-handed neutrino interacting weakly there is an electron interacting 
electrically, the two interaction processes can add up independently if the 
electron is right-handed, since right-handed electrons cannot interact with 
left-handed neutrinos. Hence a possible standard model interpretation of the 
zero 02^'' would be 

af^^ = aez(3me) + a2afc(3"^i.J = a?^^ + «ez(3m^J— ^--^ — — (52) 

4 Sm Oy/ cos^ OiY 

In the harmonic oscillator interpretation of eq. (l46l) . we choose B massless, 
/i = ul, /2 = in addition to the process already described by a^'^K 
Putting in the experimentally measured value of sin^ 9w = 0.2315, we obtain 
for the right-hand side of eq. ( l52l) the value 0.01756, which coincides to 4 
digits with the observed stationary value of the moduli field 02 = 0.01755. 

Next, let us interpret ai .In analogy to the joint appearance of u and 
e, we should also expect to find evidence for a weakly interacting w-quark, 
together with a (i-quark interacting electrically. Clearly, the w-quark could 
also interact electrically, but for symmetry reasons we expect the pair {u, d) 
to interact in a similar way as (z/, e). A right-handed u-quark interacts weakly 
with the coupling 

4 sin^ Ow 



w 



Adding up the electrical interaction strength of a (i-quark, a natural inter- 
pretation, quite similar to that of the zero 0/2 , is 

af^^ = afiiSrud) + a;^^„fc(3m„) = af^^ + ^OeiiSmu f^'^f^ (54) 

9 cos"^ Uw 

The harmonic oscillator interpretation of this moduli state is B massless, 
/i = ur, /2 = ur in addition to the process underlying a^i^\ Numeri- 
cally, taking sin^ 9w = 0.2315 and evaluating the running aei using ^ 5 
MeV, we obtain for the right-hand side of eq. (15^ 0.001800, which should be 
compared with a^^ = 0.001801. A gain we have perfect agreement within 
the first 4 digits. It is remarkable that the same universal effective value 
sin^ 6w = 0.2315 can be used consistently for both leptons (couplings a'^^\ 
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Figure 4: Interaction energy of the 2A dynamics. 



a^^^) and quarks (couplings a^'^\ af^'). Inverting the formulas, the ob- 
served stationary values of the moduli fields can be used to predict that the 
weak mixing angle is sf ^ 0.2315. 

Note that generally the backward coupling form of the = 3 chaotic 
fields seems to describe a spinless state (formed by cr and t/^, respectively 
dL and ur), whereas the forward coupling form just describes one particle 
species with non-zero spin (e or d). A similar statement will turn out to hold 
for the N — 2 theories, replacing fermions by bosons. 



,(3B)^ 



6.3 The 2 A dynamics — strong interaction strength at 
the VF-mass scale 

If electroweak coupling strengths are fixed by suitable moduli potentials gen- 
erated by chaotic fields, then something similar should also be the case for the 
strong coupling strength a^. Let us now look at chaotic fields with N = 2. 
Fig. 4 shows the interaction energy W{a) of the 2A dynamics. Only one 
stable zero is observed: 

af^^ = 0.120093(3) (55) 
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We notice that it numerically seems to coincide with the strong coupling 
constant at the W- or Z mass scale, which is experimentally measured as 
as{mz) ~ 0.118. 

For symmetry reasons, it seems plausible that if the = 3 dynamics 
fixes the electroweak couplings at the smallest fermionic mass scales, then 
the N = 2 dynamics could fix the strong couplings at the smallest bosonic 
mass scales. The lightest massive gauge boson is indeed the W^. Hence 
our physical interpretation associated with a^"^^ would he B = W^, fi = u, 
f2 = d (respectively /i = d, f2 = u), and since N = 2 formula (H6|) implies 

a'i^^ = as{mw + 'mu + md) ^ as{mw)- (56) 

Since the ly-mass is known with high precision, eq. (1561) and (1551) yield 
quite a precise prediction for the strong coupling a^. We can evolve it with 
high precision to arbitrary energy scales, using the well known perturbative 
formulas from QCD, obtaining 

a^(m^o) = 0.117804(12). (57) 

This prediction of from the zero of the chaotic 2A dynamics is clearly 
consistent with the experimentally measured value 0.118 and in fact much 
more precise than current experiments can verify. 



6.4 The 2B dynamics — the Ughtest scalar glueball 

The interaction energy of the 2B dynamics is shown in Fig. 5. W{a) has only 
one non-trivial zero 

af^) = 0.3145(1). (58) 

It is a stable zero, so it should describe an observable stable stationary state of 
a moduli field. One possibility is to interpret this as a strong coupling at the 
lightest glueball mass scale. The lightest scalar glueball has spin J^^ = 0++ 
and is denoted by gg^^^ in the following. In our oscillator interpretation we 
take B = gg^^~^, fi = u, f2 = u, thus 

a^i^^ = as{mggO++ + 2mu) ~ as{mggO++). (59) 

The 2B dynamics then describes two bosons at the same time (two gluons 
forming a glueball), similar to the 3B dynamics, which described two fermions 
at the same time (left and right handed). In both cases a spin state is 
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Figure 5: Interaction energy of the 2B dynamics. 

formed in total. In lattice gauge calculations including dynamical fermions 
the smallest scalar glueball mass is estimated as mggO++ = (1.74 ±0.07) GeV 
j21j and at this energy the running strong coupling constant is experimentally 
measured to be as ~ 0.32. This clearly is consistent with the observed value 
ot a\ ' . 

6.5 The 2A and 2B dynamics — towards a Higgs 
mass prediction 

Two chaotic field theories are still remaining, namely those with N = 2 and 
antidiffusive coupling. The interaction energies W{a) of the 2A^ and 2B ~ 
chaotic fields are shown in Fig. 6 and 7. Let us now try to find a suitable phys- 
ical interpretation for the observed smallest stable zeros a'l^ ^ = 0.1758(1) 
and af^~^ = 0.095370(1). A gain let us be guided by discrete symmetry con- 
siderations when attributing the various stationary moduli states to standard 
model particles. We saw that the smallest stable zero of the 2A dynamics 
describes a boson with non-zero spin (the lightest massive gauge boson W^) 
and the smallest stable zero of the 2B dynamics a boson without spin (the 
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Figure 7: Interaction energy of the 2B dynamics. 



22 



lightest scalar glueball) . Thus it seems reasonable to assume that the smallest 
stable zero of the 2 A ~ gives us information on yet another bosonic particle 
with non-zero spin, possibly the lightest glueball with spin J^*^ = 2"*"+, and 
the smallest zero of the 2B ~ dynamics on yet another bosonic particle with 
spin 0, possibly the lightest Higgs boson. 

Let us start with ai"^ \ Our physical interpretation of this stationary 
moduli field state in the sense of eq. (H^ is S = gg^^^, fi = q, f2 = where 
q, q are suitable quarks. From the strong coupling interpretation 

af^"^ = (60) 

the energy Ei = mgg2++ + 2m^ can again be determined from the usual QCD 
formula, evolving our previously determined as{mw) to lower energies. One 
obtains Ei ^ 10.45 GeV. In lattice gauge calculations the mass of the lightest 
2++ glueball is estimated as being roughly 2 GeV. We thus get the correct 
order of magnitude of the 2"*"+ glueball mass if we assume that the quarks 
in eq. (jlS]) are bottom quarks. With this interpretation a glueball mass 
mgg2++ = (10.45 - 2 ■ 4.23) GeV =2.0 GeV is predicted, using nib = 4.23 
GeV. 

Next, let us consider ai . Having already obtained information on 
the W boson in section 6.3, it is likely that there are also moduli fields that 
encode the Higgs boson. Our physical interpretation is B = H, fi = q, f2 = q, 
where H is the hghtest Higgs boson and q, q are suitable quarks. The strong 
coupling interpretation 

a^") = as{E2) (61) 

yields E2 = 483.4(3) GeV =mH + 2mq. However, experimental and theoret- 
ical arguments imply that the Higgs mass should be in the region 100. ..200 
GeV. Hence we only obtain a consistent value for the Higgs mass if we assume 
that the quarks involved are t quarks. This is similar to the zero a^i^ \ where 
the quarks involved were also heavy quarks. Generally, chaotic fields with 
antidiffusive couplings seem to encode information on heavy quarks rather 
than light ones. 

From the self energy of the 2B dynamics, one can obtain quite a precise 
prediction of the free top mass, namely rrit = 164.5(2) GeV, corresponding 
to a top pole mass of 174.4(3) GeV (see [2l [3j for more details). With this 
value the zero af^ ■* yields a Higgs mass prediction of 

mH = E2- 2mt = 154.4(5)GeV. (62) 
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This is a very precise prediction, the statistical error is very smalL But of 
course the main source of uncertainty is a theoretical uncertainty, namely 
whether our harmonic oscillator interpretation (in the sense of eq. P6|) ) of 

the zero ■* in terms of a Higgs boson and two top quarks is correct. For 
example, assuming that a supersymmetric extension of the standard model is 
correct, then the zero could also describe other particle states whose masses 
add up to 483.4 GeV. 

7 New physics? 

A few stable zeros of interaction energies remain that cannot be interpreted 
in the standard model context. The dynamics has yet another stable 
zero af"^^ = 0.07318. Interpreted as a running strong coupling, this would 
correspond to an energy E = |(mB + m/^ + nij^) ^ 7.85 TeV. One could 
speculate that this might describe a (rather large) supersymmetric particle 
scale. There are a few further stable zeros in the large coupling region: 
a = 0.9141 and a = 1 (3A), a = 0.3496 (35), a = 0.675 (2^") and a = 1 
{2B^). Possibly these stationary moduli states could describe gravitational 
couplings at the Planck scale. 

So far we only talked about the interaction energy of chaotic fields. But 
one can also look at the self energy as given in eq. (I39l) or fj40|) . Again a large 
number of local minima are observed that can be associated with known 
standard model interaction strengths, the energy E being again given by 
eq. (H6l) . We don't have the space to describe all the details here, but refer to 
[2], [3] . From local minima of the self energy one can get some rather precise 
predictions of fermion masses, in particular for the heavy fermions t, b, c, r. 
The 'landscape' generated by the self energy has also minima that can be 
identified with Yukawa couplings and gravitational couplings. 

The lighter fermions are much more difficult to deal with in this context. 
Since for very small couplings the self energy exhibits scaling behaviour with 
log-periodic oscillations [22], one is only able to give predictions of light 
fermion masses modulo 2, and there are also some other theoretical ambi- 
guities on how to associate the various minima with the light particles. For 
an early attempt to predict neutrino masses, see p]. Clearly more work is 
needed to eliminate the ambiguities for light fermions. 

Proceeding to higher energies (by investigating the chaotic dynamics for 
coupling constants a that coincide with running standard model coupling 
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constants couplings ai,a2,«3 at large energies), and using again the self 
energy as the relevant observable, one can also look for local minima that 
potentially could describe supersymmetric particles with masses in the TeV 
region. No such minima are found [3J. Generally the evidence for supersym- 
metric particles in the TeV region from our chaotic fields is pretty meagre, 
to say the least. One would expect to see lots of local minima in the relevant 
coupling region if the conventional ideas on supersymmetry breaking at the 
TeV scale are correct. The self energies of the chaotic fields do not provide 
any evidence for such particles: There are simply no minima in the relevant 
coupling region, the only exception being perhaps the zero a^^^^ of the 3A 
interaction energy. However, minima are found at the energy scale 10^^ GeV 
and at the Planck scale. See [5] for more details. If future accelerator exper- 
iments do not find any evidence for supersymmetric particles, a theoretical 
reason could be that they do not fit into the 'landscape' generated by the 
chaotic fields. 

Talking about new physics, new physics is of course represented by the 
existence of the chaotic fields themselves. For small couplings a, their equa- 
tion of state is close to w = — 1 hence they can account for dark energy 
in the universe. But for large couplings they can have an equation of state 
close to w ~ [4J. Thus, in principle at least, chaotic fields could account 
for both, dark energy and dark matter in the universe. 

8 Conclusion 

We have introduced chaotic scalar fields of as a model of vacuum fluctuations 
in the dark energy sector. These chaotic fields were used to generate poten- 
tials for moduli fields. We numerically observe that minima of the potentials 
lead to realistic standard model coupling strengths. The values of the fine 
structure constant, of the weak mixing angle, and of the strong coupling at 
the W mass scale are obtained with high precision and correspond to stable 
stationary values of moduli fields. Based on additional discrete symmetry 
assumptions, a value of the Higgs mass of 154 GeV is predicted from the 
chaotic field dynamics. This prediction can be experimentally tested in the 
near future. 

The theory described in this paper may be regarded as a 'tip of an iceberg'. 
It still needs to be embedded into a greater context. In particular, its possible 
relation to supersymmetry breaking mechanisms has to be clarified. However, 
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one statement can be made without any doubt: The nonhnear dynamics 
given by eq. appears to distinguish certain numerical values of coupling 
constants a that do coincide with known standard model coupling constants 
with very high precision. A random coincidence can really be excluded. In 
this way the chaotic fields can help to select the 'correct' vacuum out of an 
enormous number of possibilities, shaping the world around us in precisely 
the way we know it. 
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